HD-TVP-94-17 



Rigorous mean field model for CPA: 
Anderson model with free random variables 



Peter Neu 1 and Roland Speicher 2 * 

1 Institut fur Theoretische Physik, Universitat Heidelberg, Philosophenweg 19, 69120 Heidelberg, 

Germany 

2 Institut fur Angewandte Mathematik, Universitat Heidelberg, Im Neuenheimer Feld 294, 69120 

Heidelberg, Germany 
(September 20, 1994) 



Abstract 



A model of a randomly disordered system with site-diagonal random energy 
fluctuations is introduced. It is an extension of Wegner's re-orbital model to 
arbitrary eigenvalue distribution in the electronic level space. The new fea- 
ture is that the random energy values are not assumed to be independent at 
different sites but free. Freeness of random variables is an analogue of the 
concept of independence for non-commuting random operators. A possible 
realization is the ensemble of at different lattice-sites randomly rotated ma- 
trices. The one- and two-particle Green functions of the proposed hamiltonian 
are calculated exactly. The eigenstates are extended and the conductivity is 
nonvanishing everywhere inside the band. The long-range behaviour and the 
zero-frequency limit of the two-particle Green function are universal with re- 
spect to the eigenvalue distribution in the electronic level space. The solutions 
solve the CPA-equation for the one- and two-particle Green function of the 
corresponding Anderson model. Thus our (multi-site) model is a rigorous 
mean field model for the (single-site) CPA. We show how the Llyod model is 
included in our model and treat various kinds of noises. 
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I. INTRODUCTION 



During the last decades randomly disordered systems have gained much interest in sta- 
tistical physics. Especially, since Anderson's paper [|IJ in 1958 these systems have attracted 
many physisists due to the localization phenomenon. However, realistic multi-site models, 
like the Anderson model, are in general unsolvable beyond the one-dimensional case. Exact 
calculations are only possible in one pathological case - namely for the Lloyd model with 
Cauchy-distributed disorder. 

The Anderson model describes the hopping of an electron in a d- dimensional disordered 
lattice Z d . The hamiltonian is 

H = H + H 1 , (1.1) 
where Ho is the deterministic, translational invariant hamiltonian 

#o= E v\r-r'\\r}(r'\ , (1.2) 

r,r'£Z d 

and where the random disorder is assumed to be diagonal in the sites and independent 
between different sites, i.e. 

Hi= £ / r |r)(r| (1.3) 

r£Z d 

with f r being identically distributed, independent random variables. Although Hi is a quite 
simple operator its relation to H is complicated and in the present form the Anderson model 
is not exactly solvable. 

To circumvent this dilemma two strategies have been developed so far: One is to ap- 
proximate the multi-site model by single-site models which can be solved exactly 0-0, the 
other is to develop models which become exactly solvable in the mean field limit of infinite 
dimension d, infinite interaction range R or infinite number of angular momentum states 
n at each lattice site |8],|9|]. Among the former the most realistic ones are those single-site 
models which apply the coherent-potential approximation (CPA) ^-0; among the latter 
Wegner's n-orbital model || is most frequently studied. Wegner's generalization of the An- 
derson model consists of putting n electronic states at each site and describing the disorder 
by gaussian random matrices in the electronic states. Whereas for n = 1 this reduces to the 
Anderson model with gaussian disorder, the opposite limit n — > oo becomes exactly solvable. 
Interestingly, this solution coincides with a CPA-solution of the Anderson model where the 
single-site disorder is distributed according to Wigner's semi-circle law [^0| , pH| . 

This very fact has gained much interest in the debate about the range of validity of CPA 
and its connection with the mean field models. As mentioned by Khorunzhy and Pastur 
@ the infinite d, R and n limits do not coincide with CPA in general; however, they have 
similiar properties. 

The starting point of our investigation is the following observation: The main reason 
for the difficulty in solving the Anderson model is that the assumption of independence 
of the f r at different sites cannot be translated into a tractable relation between Hq and 
Hi. Hence, in our approximation of the Anderson model, we replace the assumption of 
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independence by a "non-commutative" independence: we assume the f r to be free. Freeness 
has been introduced in mathematics in the context of von Neumann algebras by Voiculescu 
T3] and has been extended to non-commutative probability theory by Voiculescu [|T2|-[T( 



and Speicher |L7| -[jl|. The assumption of freeness will allow us to calculate all physical 



quantities in our model exactly 

Freeness and random matrices are intimately connected with each other: Arbitrary her- 
mitian n x n-matrices randomly rotated against each other - via unitary random matrices 



- are in the limit n — > oo a possible representation of free random variables. |H,20fl. From 
this point of view our model has a mean field character and can be considered as a gener- 
alization of Wegner's n-orbital model to arbitrary eigenvalue distributions in the electronic 
level space. In contrast to this, Wegner was restricted to the distribution of symmetric 



random matrices, i.e. in the limit n — > oo to Wigner's semi-circular law [|TUyrif . It is this 



last restriction which we will show to be the reason that CPA and the hitherto considered 
mean field models coincide only for the semi-circle distribution. On the contrary, we can 
show that the solution of our model and the solution of the corresponding Anderson model 
in CPA-approximation coincide always if the disorder is distributed according to the same 
distribution in both cases. Thus our (multi-site) model is a rigorous mean field model for 
the (single-site) CPA. 

The long-range behaviour and the zero-frequency limit of the two-particle Green func- 
tion are universal with respect to the eigenvalue distribution in the electronic level space. 
Independently of the distribution of the disorder we find Wegner's result for the gaussian 
ensemble || that (i) eigenstates separated by an energy u> are correlated over a length L 
which diverges like | | x / 2 for uj — > 0, and that (ii) the two-particle Green function for ener- 
gies in opposite halves of the complex plane differing by uo approaches a constant for d > 2, 
diverges logarithmically for d = 2 and like |c<j| d / 2_1 for < d < 2. 

The paper is organized as follows: In sect. [Ij we introduce the concept of free random 
variables and outline their connection with random matrices and their description by non- 
crossing cumulants. In sect. |T| we introduce our model and calculate the one- and two 
particle Green function and the conductivity exactly. Sect. fV] is devoted to the connection 
of our model with CPA. In sect. [V] we discuss our model for various kinds of disorder, and 
show how the Llyod model is included in our model, and finally, in sect. [VTj , we summarize 
our main results. 



II. FREENESS, RANDOM MATRICES, AND NON-CROSSING CUMULANTS 
A. The concept of freeness and random matrices 



The concept of freeness was introduced by Voiculescu [[12] in order to treat non- 
commutative random variables in an analogous way as commutative (classical) random 
variables are treated by the concept of independence. From an operational point of view 
independence and freeness are nothing but rules for the calculation of mixed moments of ran- 
dom variables Xi, X 2 , if the moments of all X r are given, separately. Thus, independence 
of the X r means 



3 



(X r(1 )X r . (2 )X r(3 )...) = ( n X r{i) ) ( Y[ Xr(i)) 

i:r(i)=l i:r(i)=2 



(2.1) 



Freeness replaces this now by the following rule: 

Definition: Xi,X 2 ,... are free, if we have for all m e N and for all polynomials 
Pi(X),p 2 (X), ...,p m (X) of one variable X that 



whenever 



(pi(X r{ i))p 2 (X r{2) ) . . .Pm{X r ( m ))) = (2.2) 



(Pk(X r{k)) ) = (2.3) 



for all k — 1, . . . , m and r(k) ^ r{k + 1) for all k — 1, . . . , m — 1 (i.e. consecutive indices are 
different). 

First, one should convince oneself that this is really a rule for calculating all mixed 
moments of the X r 's. Let us consider the case of two variables X — Xi and Y = X 2 . For 
(XY) the definition yields that (XY) = if (X) = (Y) = 0. If X and Y have non-vanishing 
mean then by using ( p.2|) for the polynomials = X — (X)1 and P2(Y) = Y — (Y)l one 
easily finds (XY) = (X) (Y). Whereas this is the same result as for independent X and Y 
the calculation of (XYXY) via = ((X - (X)t)(Y - (Y)1)(X - (X)t)(Y - (Y)l)) yields 

(XYXY) = (X 2 )(Y) 2 + (X) 2 (Y 2 ) - (X) 2 (Y) 2 (2.4) 

and shows thereby that independence ((XYXY) = (X 2 )(Y 2 )) and freeness are quite dif- 
ferent concepts. Furthermore, freeness is really a non-commutative concept: If X, Y are 
free, we have (XXYY) = (X 2 )(Y 2 ) which shows, cf. (|2.4|) , that X and Y do not commute. 



Hence X and Y cannot be represented by classical c-number random variables. 

There exist a canonical representation of free random variables by special kinds of random 
matrices: Let U (n) be the ensemble of unitary n x n-matrices equipped with the canonical 
invariant Haar measure. Take two deterministic n x n-matrices A and B (e.g. diagonal 
matrices) and rotate them against each other randomly, i.e. X := A and Y := uBu* with 
u G U(n). Then, in the limit n —> oo, X and Y are free with respect to (n~ 1 tr[. . .]) av , where 
(. . .) av denotes the average over the ensemble of unitary matrices. Note that n~ 1 tr[. . .] 
gives the eigenvalue distribution of our n x n-matrices. This connection between freeness 



and unitary random matrices was first discovered by Voiculescu [14], and further developed 



by Speicher [ 20| . Another representation for free random variables with a special kind of 



distributions by deformed creation and annihilation operators will be discussed in Sect. |V|. 

Let us now check that the assumption of freeness of the f r in the hamiltonian Hi results 
in a definite relation between H Q and Hi, namely they are also free. 

Theorem 1: Let the hamiltonian H be given by (|1.1|)-(|1.3|). If the fi,f 2 ,... are free 
with respect to (. . .) ens , then H and Hi are also free with respect to (. . .). Here (...) = 
(( r o| ■ ■ ■ l r 'o))ens, independent of r , and (. . .) ens denotes the average over the disorder. 
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Proof: Consider polynomials Pi,P2, ■ ■ ■ and q±, q2, ■ ■ ■ with (pi(H )) = = (qj(Hx)) for 
all i, j. Then we have to show that 

(p 1 {H )q 1 (H 1 )p 2 (H )q 2 (H 1 ) . . .) = (2.5) 

and 

(q 1 (H 1 )p 1 (H )q 2 (H 1 )p 2 (H ) . . .) = . (2.6) 

We only treat the first case, the second is analogous. Note that (r\qj(H\)\r') = 5 r yqj(f r ). 
Then 

(p 1 (H )q 1 (H 1 )p 2 (H )q 2 (H 1 ) . . .) = ((r \p 1 (H )q 1 (H 1 )p 2 (H )q 2 (H 1 ) . . . |r » ens 
= J2 ((robi(^o)|r(l))gi(/ r(1) )(r(l)b 2 (ifo)|r(2))g 2 (/ r(2) ) . . . |r » cns 

r(l),r(2),... 

= £ (r \Pi(H Ml))(r(l)\p 2 (H )\r(2)) . . . ( qi (f ril) )q 2 (f r{2) ) . . .) ens • (2.7) 

r(l),r(2),... 

Since with H also Pi(H ) is translationally invariant, (?"obi(-^o)l r o) = (Pi(H )) = implies 
(r(i)\pi+i(Ho)\r(i)) = for all i, and we can restrict the sum to r(i)'s with r(i) ^ r(i + l) for 
all i. However for these terms we know that (qi(f r (i))(l2(fr(2)) ■ ■ -)ens = due to the freeness 
of the f/s and (qj{f r (j)))en B = (ij(Hi)) =0. □ 



B. Description of free random variables by non-crossing cumulants 

In the physics of disordered systems, usually, Green functions are calculated. This leads 
to the evaluation of mixed moments of - in our case - free random variables. The abstract 
definition of freeness ensures that all mixed moments are determined but we do not have a 
concrete formula for them, so far. An efficient machinery for concrete calculations are the 
non- crossing cumulants. 

Let Xi, X 2 , ... be free random variables. Then, we consider quantities k m (Yi, . . . , Y m ) for 
all m > 1, where the arguments Yi are non-commutative polynomials in X\,X 2 , . . .. These 
k m are called non- crossing cumulants and one way to define them is the following recurrence 
formula between the moments and the cumulants 



m— 1 

(Y 1 ...Y m )=Y, E 

p=0 

C{2,...,m} 



k p +i(Yi, Kj(!), ii(2), 



x(y 2 . . . Y i( i)_i )(%)+! . . . Fj( 2 )-i) • • • (Y i{ p) + i ...Y m ) 



(2.8) 



Starting with k\{Y\) = (Yi.), ( P-8|) may be used to determine k m {Y\, . . . , Y m ), succesively. 
The non-crossing cumulants were introduced in [B5 and further elaborated with regard to 
stochastic dynamics in ||21|| . Examples are (in an obvious notation) k 2 (l, 2) = (12) — (1)(2); 
for the special case of centered Y {(Yi) — for i = 1, . . . , 4), we have fc 4 (l, 2, 3, 4) = (1234) 
- (12)(34) - (14)(23). 



It follows from the results of [19] that also the following generalization of (|2.8| ) holds 
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(Yi . . . Y m Y m +i . . . Y m+ g) — (Yi . . . Y m ) (Y m+ i . . . Y m+ e) 



m i 



+ YI k p +q( Y i(i)i ■ ■ ■ > Y i(p)' Y j(i)> ■ ■ ■ > Y i(q)) 

p=lq=l i(l),...,i(p) j (1) • .j(<?) 

C{l,...,m} e{ra+l,-..,m+£} 

x (y . . . y (1M > (y (1)+1 . . . y (2) _ x ) . . . (y j(g)+1 . . . y m+ ,> . (2.9) 

Thus, the non-crossing cumulants give the corrections to the frequently assumed factoriza- 
tion of (y . . . y m y m+ i . . . Y m+i ) into (y . . . y m )(y m+ i . . . Y m+i ). 

To derive the connection between freeness and non-crossing cumulants we will use another 
characterization of the non-crossing cumulants (which is equivalent to ( |2.8| , [2.9| )), namely, 
they are uniquely determined by 

(Y) (2.10) 







h(Y) 


k m (Yi, ■ ■ 
i-l 

^ ] km+k- 
k=0 


■ j y%, y+i 


j • • • j Y m) 


-i(y, • • • , 


Yk, y+i, 


m—1 

e=i+i 


-£(y,..., 





Ye). (2.11) 



Note that in the first term on the rhs of ( [2.11 ) y and y +1 are multiplied with each other, 



and that in the second and third term Yfc+i, . . . , y and y+i, . . . , Ye are skipped in k m+ k-i 
and fc m+ j_^, respectively. Eq. (|2.11| ) allows to reduce all higher cumulants to k\. It is quite 
easy to derive from ( |2.10| , p7LTD the following properties: 

Remarks: 

1. k m (Yi, . . . , Y m ) is a multi-linear function in y, . . . , Y m . 

2. k m (Yi, . . . , Y m ) = for m > 2 if at least for one i we have y = 1 (of course fci(lL) = 1). 

Up to now we have not used any freeness; we have just defined non-crossing cumulants 
as special polynomials in the moments. That this definition gives us indeed the right tool 
for handling freeness shows the next proposition (for a more detailed proof see (H 



PROPOSITION 1: For each % = l,...,m let y be a polynomial in one vari- 
able X r ({) for some r(z), y = Pi(X r ^), and assume Xi,X 2 ,... to be free. Then, 
A; m (pi(X r (i)), . . . ,p m (X r (m))) — whenever there exists at least one pair i,j with i ^ j 
and r(z) ^ r(j) (i.e. such that y and Yj are free). 

PROOF: By ( |2.11 ) we can glue together neighbouring y,y + i with r(i) = r(i + 1) and 



hence we can assume that r{i) ^ r(i + 1) for all % = 1, ... ,m — 1. Next, we write again 
y = (y — (y)li) + (y)H and, by Remark 1 and 2, we can restrict ourselves to the case where 
all y are centered, i.e. (y) = for alH = 1, . . . ,m. But then we can reduce - by using 
( f2. 11| ) (or equivalently ( p.8|) ) and induction - k m (Yi, . . . , Y m ) to ki(Yi . . . Y rn ) = (y . . . Y m ), 



which vanishes by the definition of the freeness. □ 

Thus, the quite implicit definition of freeness, namely that very special mixed moments in 
free variables vanish, has now been replaced by the statement that all non-crossing cumulants 
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with at least two different free variables vanish without any restriction on (pi(X r ^)) (cf. 
the definition of freeness). So we have e.g. k 3 (Xi, Xf , X 2 ) = if Xi and X 2 are free - 
independent of the values of (Xi), (Xf) and (X 2 ) - whereas k 3 (Xi, Xf , Xi) ^ in general. 
It is this very property of free random variables which will allow us to calculate the one- 
particle and the two-particle Green function of our model in the next section exactly. Note 
that the non-crossing cumulants play exactly the same role for free random variables as the 
usual cumulants do for independent random variables. 



III. SITE-DIAGONAL ANDERSON MODEL AND FREENESS 

A. The model 

We consider now the following model of a randomly disordered system: The hamiltonian 
H is given by H = Hq + Hi , where Ho and Hi are defined in ( |1.2| ) and ( |1.3| ), respectively, 



and where we assume the f r to be identically distributed and free with respect to the average 
(. . .)ens- Due to our discussion on the relation between freeness and random matrices we 
also have the following concrete realization of the model: At each lattice site r we have n 
electronic levels \ra) numbered by a — 1, . . . , n. Let / = ((a\f\P))a,p=i be a fixed operator 
in the electronic level space and put 

f r := U r ful (3.1) 

where the u r are unitary random matrices in the electronic level space chosen independently 
for different sites r. This means that we act at each site r with a copy f r of the given 
operator /, but that the basis for f r and the basis for f T > are rotated randomly against each 
other for all pairs of different sites r ^ r'. According to our remarks around eq. ( |2.4| ), the 
f r become free in the limit n — > oo. Thus, freeness is the correct mathematical notion for 
the n — > oo-limit. 

Obviously, our model is nothing but a generalization of Wegner's n-orbital model ||. In 
his original formulation Wegner has chosen the f r = ((l/\/^)/" /3 )a,/3=i as gaussian random 
matrices at each r, such that the entries f" 13 and f^, S are independent for different sites 
r 7^ r'. Since random rotations between different sites do not alter this ensemble, it follows 
that also in this model the / r 's become free in the limit n — > oo. Hence Wegner's model is 
included special case in our model. 

One should however realize that our model is much more general than Wegner's, since the 
gaussian ensemble corresponds to choosing a matrix with Wigner's semi-circular eigenvalue 
distribution |[T0| , pT| for / in (|3.1|) . In contrast to that we are totally free in choosing any 
eigenvalue distribution for / and hence for Hi. 

It is interesting to note that freeness is already to some extend contained in the original 
Anderson model (f r being independent, random c-numbers). For instance we find in this 



case for (H HiH Hi) the same result as in (|2.4j) : 

(H HiH Hi) = ((r \H HiH Hi\r )) 

( y ] V\r— rn\frV\m— r\fi 
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^ "j V\r—ro\'U\ro~r\{fr)eris{fra)eris V(){f ro ) ens 

= ( H l) 2 (Y, V \r-r mr -r\ ~ + (#l)«0 

r 

= (Fx) 2 ^ 2 ) + {Hl){H,f - (Hi) 2 (H ) 2 . (3.2) 

Thus, the usual Anderson model yields freeness between Hq and H% for small moments, 
but something uncontrollable for higher moments which precludes the model to be exactly 
solvable. 



B. One-particle Green function 



We want to calculate the averaged one-particle Green function (IPG) defined by 

1 



G(r,r';z) :=((i 



z-H 



(3.3) 



In matrix notation this reads G(r, r';z) = n~ l Y,a(( r i a \ [ z ~ H] -1 |r', a)) ens . Let us first 
concentrate on its diagonal part G(r ,r ;z) which is independent of r due to translation 
invariance. Let us introduce the short-hand notation G(z) := G(ro,ro; z) and (...) := 

((^o| ■ ■ ■ ko))cns, thus 



G(z) 



E 

n=0 



y n+l 



t z-(H + Hi)/ : 
On calculating this quantity we assume that we know the IPG of H and Hi, separately. 

1 



(3.4) 



G (z) 
Gi(z) 



z-H 
1 

z-Hi 



~ « (r \H-J\r ) 

2.^ ~n+l ^ 



n=0 

oo 

E 

n=0 



Z" 
7 n+l 



n=0 



y n+l 



E 

n=0 



yn+l 



(3.5) 
(3.6) 



Thus, our problem consists of calculating moments of Hq + Hi given the moments of Hq 
and Hi, separately, where due to our assumption of the f r being free and by Theorem 1, H 
and Hi are free. In analogy to the usual convolution, which describes the sum of indepen- 



dent random variables, we have to calculate the so called free convolution |i3||15l1 of the free 
random variables Ho and Hi. At this point the difficulties with the usual Anderson model 
become evident: Independence of fi,f 2 ,... does not imply a definite relation between Hq 
and Hi so that no well-defined notion of convolution between the distributions of Hq and 
Hi exists. 



Theorem 2: Let the hamiltonian H be given by (|TJT])-(|0|) where fi, f 2 , 
identically distributed. Then, the diagonal part of the IPG is given by 



are free and 



G(z) = G Q z-Ri[G{z)\ 



(3.7) 
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where R\ is determined by 



The off-diagonal part of the IPG is given via the Fourier transform 

G(r,r';z) = f G(q; z) e iq{r ~ r,) (3.9) 

Jq 

by 

G(g; z) = w s 1 5 F777 , ( 3 - 10 ) 

z-v{q) - Ri[G{z)\ 

where iv_ r / = J q v(q) e ig ^~ r '\ Here, / ? := Jisz^q, V being the volume of the first 
Brillouin zone (1BZ). 

Remark: 

Note that the diagonal part G(z) entirely determines the off-diagonal part since with 
G (q; z) = [z — eqs. (|3.9|) and |3.10|) are equivalent to 

G(ry;z)=G (r,r';z-R 1 [G{z)}) . (3.11) 



PROOF: We first prove Theorem 2 for the diagonal part of the IPG. Proposition 1 implies 
that the non-crossing cumulants of H = Hq + Hi are additive, i.e. 

k rn {H) = k m {H Q )+k m {H l ) (3.12) 

with the short-hand notation k m (H) := k m (H, . . . ,H). Hence the free convolution is lin- 
earized by the non-crossing cumulants as the usual convolution is linearized by the usual 
cumulants. It remains to derive a relation between the non-crossing cumulants and the IPG. 
If we specialize ( |2.8|) to Y\ = . . . = Y m = H we obtain 

m m—p 

(H m ) = Y, E k p (H)(H j M)...(H j M) . (3.13) 

p=l j(l),...,j(p)=0 

j(l) + ...+j(p) = m-p 

If we now define 

00 

R(w) := £ w m k m+1 (H) (3.14) 

m=0 

then we find with ( |3.4| , |3.13[ ) the relation 

G < 2 > = j^km ■ (3 - 15) 
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Thus, R can be considered as the self-energy of G(z) which depends self-consistently on 
G(z) itself. Relation ( p. 15 ) and its equivalent form 



G[R{w) +W- 1 } 



w 



(3.16) 



are due to Voiculescu [13[ who calls R the i?-transform of H. The derivation given here 
using non-crossing cumulants has first been given in for a dynamical generalization see 

1 



2jfl . In the same way we can write 

G (z) 



G x {z) 



z — Rq[G q {z)) 

1 

z-R 1 [G 1 (z)] 



where 



Ro(w) := £ w m k m+1 (H ) 

m=0 

oo 



m=0 



Because of ( |3.12|) R is also additive 

R(w) = R (w) + R 1 (w) 



(3.17) 
(3.18) 

(3.19) 
(3.20) 

(3.21) 



Defining y by G{z) = G {y) we get z - R [G(z)} - Ri[G{z)} =y- R [G (y)} =y- R [G(z)} 
and thus z — R\[G(z)} = y from which we finally derive (|3.7| ) which proves together with 
( |3-1S| ) the first assertion of our theorem. Note that ( |3.7| ) reduces to ( |3.8| ) if we put H = 0: 
then Gq(z) = z~ x and G — G\. Note also that there is only an apparent asymmetry between 
Hq and Hi, since we may write in the same way 



G{z) = d z-Rq[G(z)] 



(3.22) 



Let us now treat the off-diagonal part of the IPG. Again, our assumption of freeness of 
the f r will guarantee that we can derive an exact expression for G(r,r',z). Using Dyson's 
equation gives 

G(r,r'; z) = G (r,r'; z) 

oo 

+ £ £ (G (r,ri;z)f ri G (ri,r 2 ;z)f r2 ...f rm G (r m ,r';z)) cns 

m=l ri,...,r m 

oo 

= G (r,r';z) + J2 J2 G o(r,ri; z) . . . G (r m ,r'; z)(f ri . . . f r J ens 

m=l r 1 ,...,r m 
oo 



G (r,r';z) + G (r,ri; z) . . .G (r m ,r'; z) 

m=l ri,...,r m 

TO— 1 

X k p+ i(f ri , /r i(1) , • • • > fr i[p) ) 

p=0 i(l),...,i( P ) 

C{2,...,m} 

x (fr 2 ■ ■ ■ /r i (i)_i)(/r i(1 ) +1 • • • fr i(2 )-i) • • • (fr i{p)+1 ■ ■ ■ fr m ) 



(3.23) 
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where we have used the recurrence formula (|2.8| ) for the non-crossing cumulants. Due to the 
freeness of the f r and Proposition 1 only such terms contribute where r x = = . . . = r^, 
which yields after some resummations 

G(r,r';z) = G (r,r';z) 

oo 

+ G o( r > r ii z ) k P+i(fn)G(r u r x ; z) p G(r u r'; z) 

p=i n 

= G (r, r'; z) + R x [G{z)\ ]T G (r, n; z)G(r 1 , r'; z) (3.24) 

n 

where we have used k p (f ri ) = k p (Hi) and G(z) = G(ri, rj.; z). Note that, once G(z) is 
known, ( |3.24j ) is a linear system of equations for G(r,r'; z) which can be solved by Fourier 
transformation. Then, eq. ( |3.24j ) reads 

G(q; z) = G (q; z) + R 1 [G(z)]G (q ] z)G(q; z) (3.25) 

which yields with 

Go(q;z) = 1 —, (3.26) 

z — v(q) 



the second assertion eq. (|3.10| ) of our theorem. □ 



We finally comment on the analytic structure of the solution (|3.7| ). By definition G(z), 
Gq(z) and G\(z) are holomorphic functions in the upper complex half plane C + . However, 
it is a priori not clear whether G(z) initially defined by eq. (|3.7| ) only in a neighbourhood 
of oo has an analytic continuation to C + . 

Clearly, the implicit definition of G(z) by (|3.7|) is unique except for the critical points 
z G C + where G'(z) = 0. Let us denote this set by D = {z G C + \G'(z) = 0} and by 
A = G{D) their critical values; analogously, we define Dj and Aj for j = 0, 1. Then 
Voiculescu has shown in [plj for compactly supported measures of H and Hi and G(z) 



implicitly given by fl3.7p that 

1. G(C+) C G (C+) DGi(C+). 

2. Rj[Gj(z)], (j = 0,1,2) has an analytic continuation from oo to C + (here G2 = G, 
R 2 = R). 

3. if A fl Ai = 0, the function y(z) := z — Ri[G(z)} has an analytic continuation from 
00 to C+. 

The assumption Ao fl Ai = implies that if w~ Y + Rq{w) - which is the inverse of Go (cf. 
eq. ( |3.16|) ) - has a branching point of order p > 1 at £0 then £0 is not a branching point of 
w^ 1 + Ri(w) and hence Co is a branching point of order exactly p for w~ l + R(w) |IBfl; this 



guarantees the uniqueness of the analytic continuation of eq. (|3.7|) . This assumption is not 
restrictive for practical purposes in physics as long as the measures of H Q and Hi do not 
coincide. However, in this case we know due to the freeness of H and Hi (cf. Theorem 1) 
that R(w) = 2R (w) = 2R 1 {w). 

Before we develop further our formalism for the 2PG, let us give Ri[G(z)} a clear physical 
interpretation. The diagrammatic representation of eq. ( |3.24| ) is 
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here: Diagram 1 

in obvious notation. Thus, R\\G(z)\ behaves like an effective local potential which scatters 
the propagating electron incoherently at each lattice site. The total scattered wave is the sum 
of the contributions from each lattice site without interference terms. From this viewpoint 
our model has CPA character. We will come back to this in sect. |V|. 



C. Two-particle Green function 



In this section the averaged two-particle Green function (2PG) 



G(r, s,s,r;zi, z 2 ) := ( (r 



zi-H 



8){S 



z 2 -H 



(3.27) 



will be calculated. In matrix notation this reads Q(r, s, s', r'; z±, z 2 ) = n 1 J2 a ,p(( r i a \ [ z i 
H]~ 1 \s, f3)(s', f3\[z 2 - H]~ 1 \r',a)) cns . We define the 2PG of Hi by 



Gi(z l7 z 2 ) 
From the identity 



zi - H 1 z 2 - H l 



fro z 2 fro 



(3.28) 



[zi - z 2 ) 



1 



zi-H 
Z\ - z 2 



S (S 



z 2 -H 



[ Zl -H){z 2 -H) 
= G(r,r';z 2 ) - G^r';^) 

one obtains the sum rule 

{zi - z 2 ) G{ r i s i s i r 'i z i, z 2) = G(r, r'\ z 2 ) - G(r, r'\ z x ) 

s 

which reduces for Qi to 

{zi - z 2 ) <3i{z!, z 2 ) = Gi(z 2 ) - G^zx) . 
The 2PG of our model is now given by the following 



(3.29) 



(3.30) 



(3.31) 



Theorem 3: Let the hamiltonian H be given by ( |1.1| )-( |Q| ) where fi, f 2 , ■ ■ ■ are free and 
identically distributed. Let further G(z) and G(r, r';z) be given by ( |3.7| )- (|3~T0| ). Then, for 
Im zi/2 ^ 0, 



G(r, s, s', r'; z u z 2 ) = G{r, s; z x )G(s' ', r'; z 2 ) 
+ K 1 [G(z 1 ),G(z 2 )} 

x ^ G(r, r"; zi)Q(r", s, s 1 , r"; z l7 z 2 )G(r", r'\ z 2 ) 



(3.32) 



12 



where 1Zx(wx,w 2 ) is determined by 

TZi(wi,w 2 ) = . (3.33) 

Wi - w 2 

Remark: 

Note that, once the IPG is known, ( |3.32| ) is a linear system of equations for the 2PG. 
PROOF: Again we use Dyson's equation to obtain 

g(r, s, s', r'; z x , z 2 ) = ( I G (r, s; z x ) + G o(j, n; Zx)f n . . . f rm G (r m , s; z x ) \ 

\ I m=l ri,...,r m ) 

x (Go(s',r';z 2 )+f; £ G (s', si; z 2 )f Sl . . . f Se G (s h r'; z 2 )}\ • (3.34) 
I e=is u ...,s e ) I __ a 



As before we can express this in terms of the non-crossing cumulants by using relation ( |2.9|) 
for (f ri . . . fr m f Sl ■ ■ ■ fs e )- Due to the freeness of fx, f 2 , ■ ■ ■ we can again use Proposition 
1 and therefore we can restrict the summation in ( |2.9|) to terms with r^x) = . . . = 
= = . . . = After some resummation we finally obtain 



Q{r, s, s', r'\ zx, z 2 ) = G(r, s; Zx)G(s', r'; z 2 ) 



+ E Y, k P +q (H 1 )G(r,r" ] z 1 )G(r",r" ] z 1 r 1 

p,q=l r" 

x g{r", s, s', r"; z x , z 2 )G(r", r"; z^Gir", r'; z 2 ) . (3.35) 



Defining 



Kx(wx,w 2 ) := J2 k p +q( H i) w i 



P- 1 n„<l- 1 



Wo 



p,q=l 



Yk n (Hx) £ w\ 

n=l p,<?>l 

p+q = n 



P- 1 „„<1- 1 



= Y,k n (Hx)^-^ (3.36) 
wx - w 2 

yields with ( |3T20| ) the eqs. ([T3|) and flSTBp . □ 

Following Wegner || (his eqs. (4.18)-(4.21)) one easily shows that eq. ( |3.32j ) is consistent 
with the sum rule ( |3.30| ). Using this sum rule together with ( |3.33|) yields the relation 

t> \ni \ m wS^ni \ Ri[G{zij\ - Rx[G(zx)\ , . 

Kx[G(zx),G{z 2 )\ 2^y(r,s,s,r;zx,z 2 ) = (3.37) 

s Zx — Z 2 

which closely resembles a Ward identity. Since Rx[G(zx)] is the self-energy of the IPG the 
lhs may be interpreted as a vertex function. The structure of eq. ( |3.32j ) permits still another 
characterization of lZi[G(zx),G(z 2 )] as a local effective electron-electron interaction. This 
becomes particularly pronounced if one represents (|3.32|) diagrammatically by 
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here: Diagram 2 



Thus eq. ( |3.32| ) has a ladder structure with an effective electron-electron interaction 
lZx[G(zx),G(z 2 )). This interaction is a contact interaction, i.e., after averaging, the two 
electrons propagate independently through the lattice unless they meet at same site. As for 
IPG our model resembles CPA character. We come back to this in sect. ITV. 



D. Long-range behaviour and conductivity 

Long-range behaviour. To discuss the long-range behaviour of the 2PG let us consider 
its connected part for r = r' and s — s' 

C(r, s; Zx, z 2 ) := <3{r, s, s, r; z 1 , z 2 ) - G(r, s; z 1 )G{s, r; z 2 ) (3.38) 

and its Fourier transform 

C(q; z h z 2 ) = £ C(0, r; z u z 2 ) e^ r . (3.39) 

r 

The disconnected part of the 2PG yields only a short-range contribution and will not be 
discussed in the following. 

Defining Qi 2 (r, s; Zi, z 2 ) := G(r, s; Zx)G(s, r; z 2 ) and its Fourier transform 

M<?; z h z 2 ) := £ G(r, 0; Zl )G(0, r; z 2 ) e iqr 

r 

= I G(q'; Zl )G(q'-q;z 2 ) (3.40) 

Jq' 



with G(q; z) given by ( p,10| ), and using ( |3.32j ) one obtains with the abbreviation 1Z\ = 
n 1 [G{z 1 ),G(z 2 )} 

C(q; Zx, z 2 ) = Tlx ( GMq; z u z 2 ) ) + Tlx Q\ 2 {q\ z 1} z 2 )C(q; z u z 2 ) , (3.41) 
which has the solution 

C^ ^^f;^ . ,3.42) 

Now we want to show that (7(0; zi, z 2 ) diverges if zi and z 2 approach the same energy E 
from different halves of the complex plane along the branch cut of G, that is for 

g(E) = — [ G(E - t0 + ) - G(E + i0 + ) ] ^ . (3.43) 
This follows from the decomposition 

n( \nl \ G(g; z 2 ) - G(q; Zx) 

G(q] Zl)G(q] Z2) = zx-Rx[G{zx)]-z 2 + Rx[G{z 2 )] ' (3 " 44) 
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where we have used ( |3.10|) , which yields 

MO; *.*) = ^ + Rl( „_ Jl) X\ ( G W - Gfe) ) ' (3 ' 45) 

Thus the denominator of ([3.42 ) vanishes for g = and limzi/2 = -E since the second term of 
( P.45f ) vanishes because of ( |3.43|) . From this one concludes that C(q;zi,z 2 ) has a diffusive 
pole. 

To make this explicit we take 

z 1 = E + -u , z 2 = E--u (3.46) 

where E is real and u has an imaginary part of sign s. Inserting this in (|3.45|) and assuming 
as Wegner || cubic symmetry with coordination number n one finds by expanding around 
q = in leading order for small uo and q 2 

nC(q; z u z 2 ) = f-^^r + aA (3.47) 



with 

2 



A — — 



2ttq(E) 
( fi(E)Y dG 12 (q; z u z 2 ) 



I g(E) I dq 2 



9 =o 



(3.48) 



Here we have used (|3.33|) to write 



with the measure of R\[G(z)] 

f i(E):= — [R 1 [G(E-iO + )]-R 1 [G{E + tO + )]] . (3.50) 

27TZ 

In the case of the gaussian ensemble with covariance M we will show in Sect. that 
H(E) = Mg(E), i.e. Tl x = M, so that we find Wegner's result § (his eqs. (4.40)-(4.41)). 
Thus eq. ( p.47|) merely differs from the corresponding solution for the gaussian ensemble by 
a redefinition of the constant A. The function C(q; zi, z 2 ) essentially determines the long- 
range and the uj —>■ 0-limit of the 2PG. Without further calculations one concludes from this 
that the qualitative behaviour of the long-range- and the uj — > 0-limit of the 2PG does not 
depend on the distribution of the disorder in the electronic level space. 

In more detail, following Wegner [§] (his Sect. V), one sees from (|3.47| ) with the definition 
of the wave vector 

/ — tuns \ 

«= — r-zrr- ) , Re K > 3.51 
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that eigenstates separated by an energy difference uj are correlated in phase over a length 

L = \k\- 1 = [2-kAq(E)/\oj\] 1/2 (3.52) 

which diverges as |c<j| _1 / 2 as uj — > 0. Furthermore, by Fourier back transformation and di- 
mensional analysis one finds for fixed r <C L and energies in opposite halves of the complex 
plane differing by uj that C(0, r; z\, z 2 ) approaches a finite value provided d > 2, diverges 
logarithmically as a function of uj for d = 2 and like Ici;!^ 2 ^ 1 for < d < 2, respectively, for 
uj —>■ 0. This implies that for d > the eigenstates are extended since C does not diverge 
as fast as |&>| -1 . Up to a redefinition of the constant A this is the same result as for the 
gaussian ensemble M. 

Conductivity. The 2PG determines the conductivity <j t (uj) via the Kubo-Greenwood 
relation IPl 



2e 2 f°° 

CT tM = ZV; / 

TTV J-oc 



uj- 1 



a(uj,E)dE . (3.53) 



Here, np(E) = [exp{(E — Ep)/T} + l] -1 is the Fermi distribution, Ep is the Fermi-energy, 
and a(uj,E) is the current-current or the density-density spectral function (cf. 0). 

At T = and in the dc limit uj — * the conductivity is given by the spectral function 
itself 

2e 2 

cx T=0 (u;) = —a(uj,E F ) , uj -> 0, (3.54) 

7T V 

where the spectral function a(uj, Ep) is given by the connected part of the 2PG || yielding 

'72 (9 ~ f UJ UJ 

a T ^ ( W - 0) = ^u; 2 2 — C (g; ^ + - + ^0 + , ^ - - + wO" 

2-7re 2 nA „.„ . . 
= ^^^ 2 (E F ). (3.55) 

Using the definition of A in ( |3.48| ) and the Fourier transform ( |3.10| ) this result can also be 
quoted as follows 

a T ^ o (uj^0) = ^^-^(Ep), (3.56) 
where n{Ep) is the spectral function of R\[G(z)} (cf. (|3.50| )) and 



q=0 



2 

g(E F )\ dg 12 (q; z-l, z 2 ) 



B:=A 



u(Ep) / dq 2 



(3.57) 

g=0 



The last two equations show that the dc conductivity at zero temperature is essentially given 
by the square of the spectral function of the Fourier transform of the IPG. Again this result 
differs from that for the gaussian ensemble by a mere redefinition of the constant A. From 
this one concludes that the conductivity is nonvanishing everywhere inside the band, and 
that localization cannot occur in our model. 
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IV. COHERENT-POTENTIAL APPROXIMATION 



One of the most effective approximation methods for the Anderson model is the (single- 
site) coherent-potential approximation (CPA) initially proposed by Soven || and Taylor 0. 
Its main idea can be summarized as follows One introduces an effective homogeneous 

medium with the propagator Go(r, s; z— E) with an effective potential E in which the electron 
moves and demands 

G(r,s;z) = G (r,s;z-i:{z)) . (4.1) 

In other words, CPA calculates G from an effective hamiltonian H e ff = Hq + J2 r ^\ r )( r \- 
The coherent potential E is determined in such a way that the difference between the actual 
and the effective hamiltonian, H — H e ff, produces on the average zero scattering at one site, 
i.e. the averaged single-site t matrix vanishes 

where G(z) = G(r, r; z) is given by 

G(z) = G (z - £(*)) . (4.3) 

Due to the translation invariance, t(z) is site-independent. Velicky has worked out this 
concept for the 2PG and has found the following CPA equation for the 2PG 

Q(r, s, s', r'; z 1 , z 2 ) = G(r, s; zx)G(s', r'\ z 2 ) 

+ C(z, z') J2 G(r, r"; z x )Q{r'\ s, s', r"; z u z 2 )G(r", r'; z 2 ) , (4.4) 

where G is given by flOD with E satisfying §J) and 

C( Z z>) = mt{z ' ]) (4 5) 

K ' } l + G(z)(t(z)t(z'))G{z'Y { } 

Here (t(z)t(z')) contains all contributions from repeated scattering at the same site. Thus 
the 2PG is given as the sum of single-site contributions in agreement with the general CPA 
philosophy. 

In view of ( |4.1|) and ( |4.4|) the solutions of our model (|3.7|) and ( |3.32|) have CPA character. 



This has already been realized by Wegner || and Khorunzhy and Pastur [[J for the n-orbital 
model. Wegner showed that the n —>■ oo limit of his model yields the CPA solution of the 
Anderson model provided that the f r are distributed according to the semi-circle law ||. In 
general the connection between CPA and models like Wegner's n-orbital model is not clear 
(see the discussion in 0). 

Here, we will show that the concept of freeness allows to put CPA on a firm mathematical 
basis, namely, we can prove the following theorem: 

Theorem 4: Let the hamiltonian H be given by ( |l.l|) -( |i~3D where fx, f 2 , . . . are free 
and identically distributed according to a distribution P. Then, the solution for the IPG 



17 



given by Theorem 2 and the solution for the 2PG given by Theorem 3 is identical to the 
CPA solution of the Anderson model where the fi, fz, ■ ■ ■ are independent and identically 
distributed with the same P. 

Proof: We first prove the assertion of the theorem for the IPG. We show that ( p.7|) 
and ( pT8| ) solve the CPA equations flO],fO]) if we identify E(z) = R x [G{z)). Using this to 
rewrite the CPA condition (14.21) as 



<*(*)> = E<(/r--Ri(«or> 



n— 1 



w 



(4.6) 



n=l 



where w := G(z), one finds 
<*(*)> 



w 



W 



n-1 



n=0 

oo n 



n—k 



EE " </,') (-*.(■ 

n=0k=0 YV 
oo / oo / n \ 



n— 1 



fc=0 
-r oo 

— y 

W 2 ^ 



\n=k 



W 



n—k 



W 



k-1 



k=0 (Ri(w) + 1/w 



k+1 



1 



— G 1 [R 1 (w) + l/w] 



1 

w 



(4.7) 



i.e. (t(z)) = 0. Here, we have used the equivalent form of (|3.8| ), Gi[R\(w) + l/w] = w, and 
the identity 



E 

n=k 



^n—k 



1 — a 



,k+l 



(4.8) 



To prove the assertion of the theorem for the 2PG we show that C(z, z') = 1Zi[G(z), G(z')]. 
We first calculate the average of the product of the single-site t matrices 



A(z,z') := (t(z)t(z')} 



fr 



fr ~ £' 



(4.9) 



v l-(/ P -E)Gl-(/ P -£')C?'/ ' 

where we have used the abbreviations G = G(z), G' = G(z f ), S = X(z), S' = £(z'). The 
average can be evaluated using (|4~2"|). Identifying S(z) = Ri[G(z)) one finds 

Ri[G) - R\[G'\ 



^ Z,Z ) G — G' + GG'(Ri[G'] - Ri[G]) 
which reduces with eq. (|3.33|) to 

A{z, z') = TZt[G, G'}/(1 - TZt[G, G']GG') . 



(4.10) 



(4.11) 
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Solving this equation for TZ\ yields the rhs of ( f4.5|) and therfore proves the second assertion 
of our theorem. □ 



In CPA both the IPG and the 2PG are entirely determined by the coherent potential 
E(z) which follows self-consistently from ( |4.2| ) and ( }4.3| ). In our model is given by 

Voiculescu's .R-transform R\ of the disorder. Thus the prescription ( |3.8|) for calculating R\ 
can be considered as the formal solution of the CPA equations for arbitrary disorder. 



V. DISCUSSION OF SPECIFIC DISTRIBUTIONS 

1. Deterministic noise 

Let us start with the trivial case where all f r are deterministic attaining the constant 7, 
i.e. they have a 5(f r — 7)-distribution for all r. Then, G\(z) = [z — 7] -1 , hence R\{w) = 7 
(cf. eq. ( |3~8l) ), hence 1Zi(wi,W2) = (cf. eq. (|3.33| )). This yields the following solution for 
the Green functions 

G(z)=G (z- 1 ) 

G(q;z) = (5.1) 

z-v{q) -7 

Q(r, s, s', r'\ z\, z 2 ) = G(r, s; Zx)G(s', r'\ z 2 ), 

i.e. the connected part C of the 2PG vanishes. Since for deterministic f r there is no difference 
between independence and freeness, this is also the exact solution of the original Anderson 
model. Of course it is just given by a trivial energy shift 7 of the unperturbed solution. 



2. Cauchy (Lorentz) noise: The Lloyd model 

Consider now the Lloyd model where the f r are distributed according to a Cauchy- 
distribution with parameter 7, i.e. 

dP(f r = 6) = ^de . (5.2) 

7r 7 Z + e z 

Note that moments, and thus also cumulants of Hi do not exist in this case, but nevertheless 
our main formulas for the connection between Gi(z) and Ri(w) can be justified in this case 
too [HH. We have 

G 1 (z) = - / — jr-y — ^ de = , (5.3) 

v ' 7T 7-oo 7 2 + e 2 z - e z + t'y v ' 

hence R\{w) = isj with s being the sign of the imaginary part of w. Using ( |3.33|) one sees 
that TZi(wi,W2) — if w\ and W2 are on the same halves of the complex plane, and 

2z5i T 

TZ 1 {w 1 ,w 2 ) = ^- (5.4) 

Wi - w 2 
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if u>i and w 2 are on opposite halves. Thus in the former case, we find the same result as in 
fl5.1|) with 7 replaced by 27. In the latter case the connected part C of the 2PG does not 
vanish so that one finds a finite conductivity 



0) 



2e 2 nB 
ttV 



'7 



(5.5) 



with 5 = -{dg i2 {q)/dq 2 



\q=0- 



3. Gaussian random matrix noise: The Wegner model 



Wegner's model consists in choosing the f r to be - in the limit n — > 00 - symmetric 
n x n-gaussian random matrices with the entries of f r and of f r i being independent for 
r 7^ r' . As we have explained in Sect. |1I A and III A this means nothing but that the f r 
are free. Thus Wegner's model is the special case of ours where the distribution of the f r is 
given by the eigenvalue disribution of symmetric gaussian random matrices, i.e. by Wigner's 
semi-cirle law, 



dP(fr 



1 



2-kM 



VAM - e 2 de 



(5.6) 



for e 2 < AM and zero elsewhere. The fact that gaussian random matrices are free explains 
quite naturally Wegner's observation that his model gives the same result as the CPA with 
a semi-circle distribution applied to the Anderson model. 
For the semi-circle law one has 



Gi(z) 



AM 



1 



2M 



MGi(z) 



(5.7) 



which yields Ri(w) = Mw. This means that only the second non-crossing cumulant is 
different from zero and thus we have 7^i(u>i, w 2 ) = M. This gives the following solution 



G(z) 
G(q; z) 

C(q; z 1} z 2 ) 
■ t=0 (uj -> 0) 



Gq 



MG{z) 
1 



z-v(q)-MG(z) 
M ( Gi 2 (q; z x ,z 2 ) 
1 - MQ 12 (q; z u z 2 ) 
27ie 2 nA 



V 



(5.8) 
(5.9) 

(5.10) 
(5.11) 



with A = — M 2 (dQi 2 (q) / dq 2 )\ q=0 . These formulas were found by Wegner || and later red- 
erived by Khorunzhy and Pastur II. The deformed semi-circle law (|5.8|) had also appeared 



earlier in the work of Pastur [23] as the solution for the problem of determing the eigenvalue 
distribution of a sum W + D of a symmetric gaussian random matrix W and a non-random 
diagonal matrix D. By our remarks in Sect. |I[ this latter problem is nothing but calculating 
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the free convolution of the distribution of W and of D and hence, in the light of Theorem 
1, the coincidence of Wegner's and Pastur's result appears as no surprise. 

One should also note that in the context of the free convolution the semi-circle distribu- 
tion plays the same role as the gaussian distribution for the classical convolution. This can 
be seen, for instance, from the fact that only the second non-crossing cumulant is different 
from zero for the semi-circle distribution, similiarly as only the second usual cumulant is 
non-vanishing for the gaussian distribution. For more details on the "free gaussian" and 
related topics, like free central limit theorem or free Poisson law, we refer to 0,Oj22|. 



4- q-noise: An interpolation 



In p5| we have introduced a new class of stochastic processes which interpolate continu- 



ously between classical, gaussian, random matrix, dichotomic, and Poisson processes. This 
construction can be adapted for quenched multi-site disorder as follows: Let the disorder at 
each site r be given by 

f r := a(a r + a\ ) + ^aja r (5.12) 

in terms of deformed annihilation and creation operators a r , a\ on some Hilbert space Ti. 
These operators satisfy at each site r the deformed canonical commutation relations 

a r a\ — qa\a r = \ , (5.13) 
a r |0) = 0, (5.14) 

where 1 and |0) denote the identity operator and the vacuum in 7i, respectively. At different 
sites r^r' the operators are assumed to be free, implying that fi, f 2 , . . . are free, i.e. 

a r a)/ = , r^r'. (5.15) 

The deformation parameter q is real and varies continuously in the interval —1 < q < 1. 
For £ = 0, the limiting cases q = 1 and q — — 1 describe gaussian and dichotomic disorder, 
respectively, whereas the case q = corresponds to Wegner's n-orbital model. The £ a\a r - 
term allows to include Poisson-like disorder [23 , 25 ] 



Eq. (|5.15|) is an alternative representation of those free random variables which can 
be represented by deformed creation and annihilation operators. To make the construction 
clear, let us formulate the original Anderson model for gaussian site-diagonal disorder (q = 1, 
£ = 0) in this language: f r := a(a r + a\) with a r a\, — a\,a r = S r yl for all r, r', in particular, 
d r o} r i — o} r /(X T for r^r', which is clearly different from ( |5.15| ). 

We can now identify moments of the random variables f r in ( |1.3|) with the Hilbert space 
vacuum expectation values of products of f r by means of a generalized Wick Theorem p5| , |27| . 
By using the partial cumulants we have calculated the IPG of Hi (cf. eq. (36) in ||25|| ) 



Ri[G!(z)] 
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1 



(T 2 g (0) 



z - £q 

z - iq 



z - - 



m ^ (3) 



l/(z - (a 2 q {0) /(z - £g (0) - (cr 2 q^/(z - £g (1) 

- (o*q®/(z - £g< 2 > - (o*q®/ ...)...) (5.16) 



where 



1 - n k+1 

q (k) . = i + q + q 2 + _^ + q k = 1 (517) 

1 - q 

Thus, 

Ri[G{z)] = G{z)o*qV>/(l + G(z)(Ri [G{z)\ - ft®) 
+ (G(z)Vg«/(l + G{z){R 1 [G{z)\ - £g«) 
+ (G(z)VgM/(l + G(z) {Rx[G{z)\ - ^ (2) ) 

+ (G(,)W 3 V •••)•••) ( 5 - 18 ) 



which together with 



G(z) = Go *-.Ri [£?(*)] (5.19) 



is a closed set of nonlinear self-consistent equations for the IPG and Ri(w). The continued 
fraction (|5.18|) can be summed in closed form for q — — 1, where one finds, for £ = 0, 



i*i[G(*)] = (^/l + Aa 2 G 2 (z)-l)/2G(z), and for g = 0, where one finds, for £ = 0, #i[G(z)] = 
a 2 G(z) and thereby Wegner's model with a 2 = M. 

VI. SUMMARY 

In this paper we have applied the concept of free random variables, invented by Voiculescu 
in a mathematical context, to the tight-binding hamiltonian with site-diagonal disorder of 
an electron in a periodic solid. The difference of our model to the usual Anderson model 
lies in the fact that instead of assuming the disorder to be independent at different lattice 
sites, we have assumed it to be free. 

Both freeness and independence can be considered as a rule for calculating mixed mo- 
ments of random variables. In contrast to the case of independent disorder, free noise does 
imply a treatable relation between H and Hi. they are also free. This finally allows to close 
the infinite hierachy of equation of motion and to calculate all physically relevant quantities. 

In sect. [I Il| we have argued that in the limit n — > oo free noise can be represented by 
n x n-random matrices which are randomly rotated against each other at different lattice- 
sites, starting from a matrix with fixed but arbitrary eigenvalue distribution in the limit 
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n — > oo. In this sense, our model is an extension of Wegner's n-orbital model for the 
gaussian ensemble to arbitrary eigenvalue distribution in the energy level space. One of the 
most striking property of our model is that both the long-range behaviour and the zero- 
frequency limit of the 2PG are universal with respect to the eigenvalue distribution in the 
energy level space. 

One surprising feature of the Wegner model is that its solution coincides with a special 
CPA solution. This generalizes also to our model. In sect. [TV] we have shown that our 
solution for the IPG and the 2PG also solves the CPA equations for the Anderson model 
with the same distribution of disorder. Note that we specify our model rigorously in the 
beginning and that we are able to calculate all quantities without any further approximation. 
Thus our multi-site model is a rigorous mean field model for the usual single-site CPA. The R- 
transform of Voiculescu, R±, may be considered as the formal solution of the CPA equations 
for arbitrary disorder. It posseses the physical interpretation as an effective local potential; 
the corresponding quantity for the 2PG, IZi, can be considered as an effective local electron- 
electron interaction. Cleary as seen in sect. [V|, aside from some specific distributions, both 
functions cannot be calculated analytically. However, due to the mentioned universality, 
the Wegner model is exemplary within many respects and it might be sufficient to restrict 
oneself to this case in the general frame of a mean field approximation. 

Furthermore, our description using the theory of free random variables and the notion 
of non-crossing cumulants allows a straightforward generalization to the case of dynamical 
disorder and thus promises to give a rigorous model for dynamical CPA. These subjects will 
be pursued further in forthcoming investigations. 
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